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ABSTRACT 

High-precision radial- velocity techniques, which enabled the detection of ex- 
trasolar planets are now sensitive to relativistic effects in the data of spectroscopic 
binary stars (SBs). We show how these effects can be used to derive the abso- 
lute masses of the components of eclipsing single-lined SBs and double-lined SBs 
from Doppler measurements alone. High-precision stellar spectroscopy can thus 
substantially increase the number of measured stellar masses, thereby improving 
the mass-radius and mass-luminosity calibrations. 

Subject headings: binaries: close — binaries: spectroscopic — celestial mechan- 
ics — methods: data analysis — relativity — techniques: radial velocities 



1. Introduction 



Precise radial-velocity (RV) measurements, with long-term precisions of a few meters 
per second, are now routinely obtained by several telescopes around the world. The most 
notable scientific achieve ment of precise RV measurements has been the detection of planets 
orbiting solar- type stars (IMayor fc Quelozl Il995l : iMarcy fc Butlerl Il996l ). In this Letter we 
suggest another application of high-precision RVs, namely, the d etection of relativistic effects 
in the Doppler shifts of close spectroscopic binary stars (SBs). iKopeikin &: Ozernoyl (119991 ) 
have already detailed the relativistic effects one expects to find in the Doppler measurements 
of binary stars. Here we focus on the effects that we expect to measure in SBs, and identify 
the additional information they provide. 

The typical velocities of components of close binary stars can be as high as 150 km s -1 , 
(3 = v/c ~ (9(10~ 4 ). The classical Doppler shift formula predicts a relative wavelength shift 
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AA/A of order f3. The next order corrections are of order f3 2 ~ O(l0~ 8 ), which translate to 
0(lms _1 ). Terms of order /3 3 are beyond foreseen technical capabilities. We thus limit our 
analysis to 0{(3 2 ) effects — the transverse Doppler shift (time dilation) and the gravitational 
redshift. In principle, long term monitoring may reveal higher-order secular terms, such as 
the relativistic periastron shift or period decay through gravitational- wave radiation, but 
here we focus only on the periodic effects. These can be detected during relatively short 
observing runs, assuming the RV measurements are precise enough. 



Recently, IZucker et al.l (120061 ) have shown that the very same effects should be detectable 
in the stellar orbits around the black hole in the Galactic Center, after a decade of observa- 
tions. The context we examine here is different, and the information that can be extracted 
may contribute to the statistics of close binary orbits. 



2. Single-lined spectroscopic binary 



The Keplerian RV curve of a single-lined spectroscopic binary (SB1) can be presented 



as: 



Vr\ = K\ cos to cos v — K\ sin to sin v + eK\ cos ui + Vi 



HO 



where K\ is the primary star RV semi-amplitude, u is the argument of periastron, e is the 
eccentricity, Vr is the center-of-mass RV, and v is the time-dependent true anomaly. The 
customary procedure to solve an SB1 is to fit this orbital model to the observed RV data. This 
fit is achieved through some optimization algorithm that scans the (P, T, e) space (period, 
periastron time and eccentricity). For each trial set of values for these three parameters the 
algorithm produces the corresponding v{t) and then solves analytically for Kqi{— K\ coscu), 
Ksi (= —Ky sinu;), and Vrq, which appear linearly in the expression for Vri- 

In order to incorporate relativity into the observed RV cu rve of an SB1 we can use 
the models developed for analyzing binary pulsar timing data. iTaylor fc Weisbergl (119891 ) 
present a detailed timing mode l of a relativistic binary p ulsar, based on the relativistic 
celestial mechanics developed by lDamour fc Deruelld (119861 ). Besides the transverse Doppler 
shift and the gravitational redshift, relativity also introduces the Shapiro delay, periastron 
advance, and period decay t h rough gravitational wave radiation. The formulae for the pulse 
delay in ITaylor fc Weisbergl (119891 ) can be transformed to the RV domain by taking their 
time derivative. This calculation shows that the only terms of order f3 2 are the gravitational 
redshift and the transverse Doppler shift, corresponding to the so-called 'Einstein delay' in 
the pulsar timing model. 



We now derive the two terms in a more didactic, albeit less rigorous fashion. In the 
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center-of-mass frame, we can use energy conservation in the classic Keplerian solution to 
relate the transverse Doppler term to the radius vector of the observed component, r\: 

V = - Gm * (L - ±) (2) 

where a± is the orbital semi-major axis of the primary orbit, and c is the speed of light. 
After some algebraic manipulation, we can estimate the corresponding modification to the 
measured RV from the transverse Doppler effect: 

K 2 ( 1 - e 2 \ 
AV TD = — V 1 + ecosz/ — (3) 

where % is the orbital inclination. 

The gravitational redshift caused by the potential of the secondary component is in- 
versely proportional to the separation of the two components: 

AA Gm 2 



A c 2 (r 1 + r 2 ) 

and the corresponding RV modification term from the gravitational redshift is: 



(4) 



AV GR = Kl {K ] t K2) (1 + eoosi/) (5) 
csm i 

After transformation to the observer frame, an additional term appears, related to the center- 
of-mass motion 

-^-K 1 (cos(v + u) + ecosu) + — , (6) 
c 2c 

where V is the magnitude of the full center-of-mass velocity vector. 



In total, relativity adds the following 0(f3 2 ) term to the measured Vri: 
K x 1 



AVm = ; 2~7 

c sin % 



1-e 2 



e{2K x + K 2 ) cos v + (2K, + K 2 ) —K x 



+ (7) 



-Ki (costucosz/ — sin a; sin i/ + ecoscjj H 



c v '2c 
The resulting expression for the modified measured RV, VJn, can now be easily simplified to: 

= K'ci cos v + K'si sin ^ + eK ci + ^ro , (8) 

by collecting together the constant terms, the terms proportional to cos u, and those pro- 
portional to sin v, and introducing the modified linear elements: 

K' S1 = -Ki + — J sino; (9a) 
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K'ci 



1 + 



V R0 



cos a; + 



e 2K! + ir 2 " 
sin 2 i c 



sin i c \2 J 2c 



(9b) 
(9c) 



Equations [T] and [8] share exactly the same structure, and thus we can still apply the 
same fit procedure. However, the linear elements are more difficult to interpret now. The 
three quantities K' S1 , K' C1 , and V R0 depend on the six elements K\, K 2 , u>, smi, Vro, and Vq. 
Thus, Equations M are under-determined and we cannot completely infer the six elements 
above, unless some additional independent information is available, or further assumptions 
are introduced. 

Such independent information may be available through precise photometry of eccentric 
eclipsing binaries. There, the shapes and widths of the eclipses as well as the phase differences 
between primary and secondary eclipses can be used to estimate ui and sini. In this case, 
we may derive K 2 - the RV amplitude of the secondary: 

2K' S1 sm 2 i f K' \ 

K 2 = — — -c I cos lu -\- —r^p- sin u ] . 10 

sin a; e \ K' S1 J 

In the above equation we neglected Vrq, as it contributes only higher order terms. By 
obtaining K 2 we effectively turn the binary into a double- lined spectroscopic binary (SB2), 
in which both K\ and K 2 are measured. Together with the known inclination, we then obtain 
full knowledge of the component masses. 

Curiously, another result of including the relativistic terms is that an eccentric binary 
should always display an apparent RV signature, even in the extreme case where the inclina- 
tion is exactly zero and the orbit is observed face on. Since Kj = 2iraj smi/ (Pa/1 — e 2 ), sini 
will cancel out in Equation [9bl Then K' CI will be finite, while K' SI = 0, and the apparent 
argument of periastron will be exactly 0° or 180°. Thus, all the RV planet candidates whose 
arguments of periastra are close to these values, can in principle be binary stars observed 
exactly face on. Nevertheless, such small values of the inclination are extremely rare and 
this possibility is not realistic. 

Equation [9c] does not contribute any new useful information, since systematic effects 
in the measurement process such as spectral template mismatch are probably larger than 
the relativistic effects. In addition, the spectra are subject to gravitational redshift by the 
potential of the emitting star itself, and the typical uncertainties regarding its mass and 
radius are also larger than the effects we discuss here. 
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3. Double-lined spectroscopic binary 



In the case of a Keplerian SB2, there are two sets of measured RVs. The two sets of RVs 
share the same fundamental orbital elements (P, T, e, u, and Vrq), and the only difference 
between them is their amplitudes K\ and —K 2 . The common procedure is to scan the space 
of the four parameters (P, T, e, us) and then solve analytically for the three linear elements 
Ki, K 2 , and Vrq. However, when we incorporate the relativistic corrections above, we see 
that we now have two RV curves with two different sets of derived amplitudes, arguments 
of periastra, and center-of-mass velocities. The two RV curves still share the same period, 
periastron time and eccentricity. We then have the following equations, corresponding to 
Equations EH and [9b] 



Vro appears in the four equations always divided by the speed of light. Thus we can 
safely use its approximate derived value, from either set of measured velocities , since the 
discrepancy will be only 0(j3 3 ). We are left with four equations with four unknowns: K\, 
K 2 , u), and sini Solution of this set of equations will yield a more accurate estimate of the 
first three unknowns, but more importantly, it will yield an estimate of sini. Retaining only 
leading order terms, we can arrive at the following solution: 




(lid) 



(lib) 




(11c) 



3e K'+K', 



si 



(12) 



uj 2 — uj[ c 



where 




(13) 



Thus, we see that relativity causes the apparent arguments of periastra to differ. Note 
that classically, there is no way to estimate sin i from pure RV data alone. The value of sin i 
is usually obtained only from the analysis of SBs that are eclipsing or astrometric binaries. 
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Discussion 



We present here an approach to extract more information from RV data of a spectro- 
scopic binary, based on the inclusion of relativistic effects. The practical potential lies in the 
time-dependent parts of the relativistic terms which are closely linked to the variation of the 
distance between the binary components. Thus, these terms are especially useful for orbits 
that are eccentric enough. 

The approach is mainly useful in the case of SB2s, but precise photometry or astrom- 
etry may add the required informa tion to utilize relat ivity in S Bls as well. Precise pho- 
tometry space missions like M OST jWalker et alil2003l ). CoRoT jBagfinl booJ ) . and Kepler 
( IBasri. Borucki &: Kochll2005l ) may be able to provide precise enough measurements of uj and 
the inclination of eclipsing binaries from the analysis of their light curves. 

In any case, it is essential that the data quality be high enough to be sensitive to 
variations of the required order, namely one meter per second or less. Currently, the best 
precision is obtained by the ESO HARPS fiber-fed echelle spectr ograph, where the RV error 
can be as low as 1ms -1 in certain cases and maybe even less (ILovis et al.ll2005l ). In the 
future, much better pr ecisions can be hoped for, on instruments designed for the "extremely 
large telescopes" (e.g.. iPasquini et al.ll2006l ). 



To demonstrate the relevance of the suggested approach in real-life cases we chose to 
examine the SB2 12 Boo. Recently, iTomkin & Fekell (120061 ) published a precise solution 



of 12 Boo based on RVs obtained at the 2.1-m telescope at the McDonald Observatory 
and at the Coude feed telescope at Kitt Peak National Observatory, with RV precisions 
of 0.1-0.2 km s _1 . The system has a period of 9.6 days, eccentricity 0.2, and both RV 
semi- amplitudes are close to 70kms _1 . These orbital parameters translate to an expected 
relativistic amplitude variation of about 10 ms" 1 (Eq. [7j). Furthermore, with a declination 
of +25°, the star is observable by HARPS. Its brightness (5th magnitude) and spectral type 
(F9IV) make it fairly r easonable to exp e ct a p recision of lms" 1 with HARPS. We used the 
available 24 RVs from ITomkin fc Fekell (120061 ) and augmented them with only 3 simulated 
HARPS measurements (assuming errors of lms" 1 ) , including the relativistic effects. For 
each assumed value of sini we produced 1000 sets of simulated measurements and solved for 
the orbital elements, using Equation [T2l to estimate sinz. Figure [1] shows the median of the 
derived values in solid line, and the 25% and 75% percentiles in dashed lines. The figure 
demonstrates that with reasonable efforts, sin i can be measured satisfactorily. In the worst 
case where sinz = 1, the standard deviation of the derived inclination is 0.14 and a few more 
precise measurements can reduce this value significantly. An additional advantage of this 
test case is that the inclination of 12 Bo o has already been measured by interferometry and 
is known to be 108° ( jBoden et al.ll2005l ). Thus if this test is performed, the derived sinz can 
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be compared to the known value. 



In real observations, more than three precise measurements may be needed in order to 
account for differences in zero points between instruments. Furthermore, The results depend 
crucially on the precision of u[ and u>' 2 , and a relative error of e in uj' 2 — u>[ will be translated 
to a relative error of 1.5e in the absolute masses. The few p recise meas urements can be 
scheduled to optimize the precision of those two elements (e.g. iFord II2006I ). 



Currently, precisions of a few meters per second are still difficult to obtain and besides 
using the best instruments available, there are also several limitations imposed by the star 
itself. Thus, early-type stars or rapid rotators, where the spectral lines are significantly 
broadened, do not lend themselves easily to high-precision RV measurements. Stellar os- 
cillations and star spots are also a concern as they can cause apparent RV modulation. In 
addition, analyzing SB2s with the same level of precision as SBls has not been easy until 
recently, when TODCOR JZucker fc Mazehl|l994f) was applied successfully to high-precision 
spectra by several teams (jZucker et al.l 12004 ; lUdry et al.ll2004t iKonackil 120051 ) . 



The effects we have examined are most useful when the orbits are eccentric and the 
RV amplitudes are large enough. Large RV amplitudes are usually typical to close binaries, 
which are expected to have undergone orbital circularization and usually have vanishing ec- 
centricities. However, eccentricity somewhat increases the RV amplitude, and even relatively 
wide binaries, with high enough eccentricities, can display quite large RV amplitudes. 

Care must be taken to model correctly any other effects of order f3 2 that might con- 
taminate the data. One such effect is the light-travel-time effect. This effect can be easily 
approximated to the relevant order by adding the following term to Vri (and a corresponding 
one to Vr 2 ) : 



AW T = — 



sin 



[v + to) (1 + e cos v) 



(14) 



An effect which should be analyzed carefully is the tidal distortion of the stellar com- 
ponents, in particular close to periastron. This distortion may affect the spectral lines, 
introducing line asymmetry, which can bias the estimated Dopple r shift. RV extrasolar 



planet surveys use the line-bisector analysis (e.g.. iQueloz et al.ll200ll ) to quantify such time- 
dependent asymmetries. Further development of this technique may be the key to disentangle 
the tidal distortion and the relativistic effects. 

One important application of the proposed method is to calibrate the low-mass end of 
the mass-luminosity relation, to better understand the stellar-substellar borderline. This 
mass regime is still poorly constrained, since low mass SB 2s are quite rare due to the special 
photometric, spectroscopic and geometric requirements (jRibasI 120061 ) . Large efforts are in 
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progress to obtain accurate stellar masses in this regim e, including adapt ive optics, inter- 
ferometry and in the future space interferometry (e.g., Henry et aliboOol ). We propose a 



new, relatively accessible tool to accomplish this goal, where the only requirements a r e spec - 



troscopic. Precise RVs for low- mass SB2s were already measured by iDelfosse et al.l (119991 ). 
Using the method presented here, their absolute masses may be derived with a relatively 
small observational effort. No other method exists yet to derive this information purely from 
RV measurements. 
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input sin / 

Fig. 1. — The results of simulations of 12 Boo including relativistic effects, assuming RV 
precision of 1ms -1 . The plot shows percentiles of the derived shu against the input simu- 
lated value. The dashed lines represent the first and third quartile, and the solid line is the 
median. 



